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Abstract. We extend results of Dos Santos Ferreira-Kenig-Sjostrand- 
Uhlmann ( arXiv:math.AP/0601466 I to less smooth coefficients, and we 
show that measurements on part of the boundary for the magnetic 
Schrodinger operator determine uniquely the magnetic field related to a 
Holder continuous potential. We give a similar result for determining a 
convection term. The proofs involve Carleman estimates, a smoothing 
procedure, and an extension of the Nakamura-Uhlmann pseudodifferen- 
tial conjugation method to logarithmic Carleman weights. 



1. Introduction 

Let C R*^, n > 3, be a bounded simply connected C°° domain with 
connected boundary. We will consider two inverse problems in fi. The first 
is the inverse problem for the magnetic Schrodinger operator, defined by 



where Dj = l^,W e L°^(0; C"), and V G L"(17). Provided that is not 
a Dirichlet eigenvalue of Hwy in ^) we define the Dirichlet-to-Neumann 
map (DN map) formally as the magnetic normal derivative 

Nw,v : / ^ (V + iW)u ■ i^lan, 

where u is the outer unit normal to and u £ H^{Q) solves the equation 
Hwyu = in O with u\qq = f. Using a weak formulation, Nwy is well 
defined bounded map from H^/^{dn) into H-^/^{dn), see also Section El 
below. 

We consider Nwy as the boundary measurements related to the op- 
erator Hwy, and the inverse problem is to determine the coefficients of 
Hwy from partial knowledge of N^y. There is gauge equivalence: one has 
N]y+\?py = Nwy when p £ W^'°°{^1) and p\qq = 0. Thus one may only 
expect to recover the magnetic field dW and the electric potential V from 
the boundary measurements. Here dW is the 2-form d{Y^^^-^ Wj{x) dxj). 

Let xo G R" \ ch(r2), where ch(r2) is the convex hull of 0,. We define the 
front face of dQ relative to xq by 

F{xo) ={xedn: {x- xo) ■ iy{x) < 0}, 

and we take F to be an open neighborhood of F{xo) in dil. Also, let 
Wtan = W — (W • be the tangential component on dO, of a vector field 
W. 

I 
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Our main result states that the magnetic field and electric potential are 
determined by measurements Nyyy made on the particular subset F. 

Theorem 1.1. Let Wj £ C^(lI;C"), e > 0, and Vj £ L'^in) for j = 1,2. 
Also assume that is not a Dirichlet eigenvalue of Hwj,Vj in ^- If 

Nw^,vJ\p = Nw,,vJ\p for all / G H^'^d^), 

then dWi = dW2 in n and (T4^i)tan = (W2)tan on dQ. If further Vj G 
then Vi = V2 in ^■ 

This theorem was proved for W G C'^{Q) and V G L°°{Q) in jS], following 
the W = case in fPJ (see also 0)- These articles introduce important 
ideas, and the main tools are Carleman estimates, the construction of special 
solutions to the Schrodinger equation, and analytic microlocal analysis. The 
special solutions generalize the exponentially growing solutions in |2Uj . 

We will employ the methods of ^ and a smoothing argument from |15| 
to work with W G W^'"'{n)nC%Ti). The assumption W G ensures 
that the Carleman estimate construction gives special solutions in H^{Q). 
To remove this assumption, we will improve the construction of solutions 
by combining Carleman estimates with the pseudodifferential conjugation 
method from jl3j . The conjugation method was given in for linear 
Carleman weights, and we extend this method to logarithmic Carleman 
weights. This is a microlocal argument using some ideas from jB] and 
and it gives the theorem for W G C^(r2). We also allow complex coefficients 
and give a more precise proof when 17 has complicated geometry. 

Earlier work on the magnetic inverse problem for n > 3 has dealt with the 
case where the full DN map is known. The uniqueness question was solved 
in for small W and in jl2j for smooth W without smallness assumptions. 
The regularity of W was relaxed to in |21j and to Dini continuous in jl8j . 
A constructive procedure for recovering dW and V from N^y is given in 
jl7j . The related inverse scattering problem has been studied in j7j . We also 
mention ^ and which consider boundary determination and stability 
for magnetic Schrodinger operators, and ^0] which proves a partial data 
result for the nonsmooth conductivity equation. 

The other problem we consider is to determine a convection term from 
boundary measurements. If G L°°(r2;R") is a convection term and if 
u G H^{0,) solves (-A + 2W ■ V)n = in with u\dn = f, we define the 
related DN map formally as the normal derivative 

Nw ■■ f ^ duu\an- 

We remark that since W is real valued, is not a Dirichlet eigenvalue of 
—A + 21^ • V by the maximum principle. The map Ny\r is again bounded 
from H^/'^{dn) into H-'^/'^{dn). The inverse problem of determining W 
from Nw can be reduced to the case of the magnetic Schrodinger equation, 
and we obtain the following result as a corollary to Theorem ll.il 

Theorem 1.2. Let Wj G (7^(11; R"), e > 0, and V-Wj G L°°{n) for j = 1, 2. 
If 

NwJ\f = NwJ\f ioT all f £ H^^dn), 
then Wi = W2 in n. 
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Earlier work on the convection term problem for n > 3 includes ^ which 
solves the full data problem for smooth W, and ^S] which gives the same 
result for Lipschitz continuous W. 

The paper is organized as follows. In Section 2 we fix some notation and 
preliminaries. Section 3 contains the facts on first order elliptic equations 
which will be needed for the pseudodifferential conjugation method. The 
Carleman estimates are in Section 4, and in Section 5 we establish the pseu- 
dodifferential conjugation method for logarithmic Carleman weights. Sec- 
tion 6 gives the construction of special solutions to the magnetic Schrodinger 
equation, and Theorems 11.11 and 1 1 . 21 are proved in Sections 7 and 8, respec- 
tively. 

2. Preliminaries 

First we fix some notation. We will always take $7 to be a bounded simply 
connected domain in R", n > 3, with connected C°° boundary. Recall that 

is simply connected if it is connected and has trivial fundamental group 
(hence, the first de Rham cohomology vanishes). Let u be the outer unit 
normal of $7, and d^u = Vu • the normal derivative of u on dQ. 

We write C(J7), C^(0), and C°°{Q) for the sets of complex valued func- 
tions which are continuous, e-Holder continuous, and infinitely differentiable, 
respectively, in Q. The corresponding spaces of compactly supported func- 
tions in n are Cc(0), C^{n), and C^{n). The notation C(0; C") denotes 
the space of n-dimensional vector fields whose components are in C(0), sim- 
ilarly for C^(r2;C"'), LP(Q;C"') etc. For real-valued vector fields we write 
C'{n;W) etc. 

li u,v G we write {u\v) = J^uvdx, and also {u\v) = J^u- vdx 

when u,v E L"^ {^1,; C^) . The L^(0) norm is denoted by ||n|| = {u\u)^^'^. If 
f,g G L^((917) let {f\g)dn = fgnfgdx. We will use the same notations for 
the pairing of distributions and elements in the dual space. 

Denote by H'^{Q) the Sobolev spaces in Q, and by W^''^{i^) the usual 
Sobolev spaces. We will need the following well-known results concerning 
multiplication in Sobolev spaces. 

Lemma 2.1. If a G L^{Q), n> 3, then au maps H^{Q) to and 

l|a''^llL2(n) < C'll«llL"(n)lhll/fi(f7)- (1) 

If furthermore a G VF^'" n L°°(J7) then u ^ au \s bounded on H^{yt) for 
— 1 < s < 1, and we have 

I|a^ll_ff=(f7) < C'l|a||vyi>"nL°°(n)lhlli?''(n)- (2) 

Proof. The first statement follows from the Holder inequality and Sobolev 
embedding. The Leibniz rule and (pQ) give for s = 1, the case s = —1 
follows by duality, and the full result is obtained by interpolation. □ 

Next we note that when W is Holder continuous, one may always assume 
V ■ W = Q after a gauge transformation which does not alter Nw,v or the 
Holder continuity of W . 
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Lemma 2.2. li W G C=(Q;C'"), then there exists p £ C^+%n) such that 
V -{W + Vp) = and p\on = 0. 

Proof. It is enough to choose p to be the solution of Ap = —V • W with 
p\dn = 0, and this Dirichlet problem has a C^'^'^ solution by jSl Section 
8.11]. □ 

We proceed to give a precise definition of N^y- Let 

H^i^) = {u£ H^{n) : Au G L'^{n)} 

with norm = + ll^^llL2(r2)- This is a Banach space 

and C°°(r2) is a dense set. The trace operator is bounded on H^{Q) with 
values in H^/'^{d^), and also u i— > d^u is bounded on H^{Q) with values 
in H~^/'^{di}). This follows by writing u G H/^{^}) as u = uq + ui where 
Alio = with uolgQ = u\Qfi, and Aui = An with uiI^q = 0. If V • = 
then for / G H^/'^{dfl), the solution u to Hwyu = with boundary value 
/ is in H/\{Q). Therefore, the DN map Nwy : / ^ (V + iW)u ■ iy\dn is 
bounded from H^/'^{dVt) to H-^/'^{dn). 

Finally, we will need the Green identity in the following form. 

Lemma 2.3. Let W £ C{Tl; C") with V ■ W = 0, and let V G U'{n). If 
u G Ha{^) n H^i^) and v G Ha{^), then 

{Hwyu\v) - {u\H^, yv) = -{d^u\v)dQ. 

Proof. Let Uj,Vj G C°°{0,) with ^ u, vj — > t; in H/^{^}). Then 
{HwyUj\vj)-{uj\H^,yVj) = {uj\d„Vj)Qn-{d^Uj\vj)3n-2i{{W-iy)uj\vj)gn. 
The claim follows by taking limits. □ 

3. Elliptic equations of first order 

To extend the Nakamura-Uhlmann pseudo differential conjugation method 
to logarithmic Carleman weights, we will need to solve first order elliptic 
equations with variable coefficients in domains in T*R"'. This is because the 
symbols of the conjugating operators arise as solutions to such equations. 
The main result in this section is ProDosition l3.41 which shows the solvability 
of a first order equation related to the logarithmic Carleman weight. 

The standard reference for the following facts is [H]. See also jH] for the 
specific case of limiting Carleman weights, which will be the case of interest 
for us. Thus, let ip = (p{x) be a smooth real function in an open set V C R" 
with V(p 7^ in y, and write P = {—h? A)e~^ . The semiclassical Weyl 
symbol of P is p = a + i6, where 

a{x,i)=e-{'^v)\ h{x,i)=2V^-i. (3) 

We say that (/? is a limiting Carleman weight (for the Laplacian) if the 
Poisson bracket {a, 6} = V^a • V^^ft — V^a • satisfies 

{a, 6} = when a = 6 = 0. (4) 

This implies that {a, h} = ca + db for some smooth c, d, and then on the set 
S = {a = 6 = 0} C r*y the Hamilton vector fields satisfy [Ha, Ht] = cHa + 
dH^. Recall that the Hamilton vector field of / is Hj = V^f ■ Vx — Vxf ■ Vg. 
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It follows that S is an involutive manifold of codimension 2, and the 
Frobenius theorem states that at each point of S there are local coordinates 
in which Hp = Ha + iH), becomes an elliptic operator + ^^2(2/)^- 

This implies local solvability of HpU = f. 

For our purposes local solvability is not enough, since we will need to 
solve a related equation in a full neighborhood of S. To do this we follow jH] 
and find a smooth function m satisfying {ma, mb} = in a neighborhood 
of E, so that Hma and H^b commute near S and in some new coordinates 
Hmp becomes ^ + 

The first lemma states that in the present case where a and b come from 
a limiting Carleman weight, there is an explicit choice for m. 

Lemma 3.1. If m = \\/ (p{x)\~'^ , then {ma,mb} = in T*V. 

Proof. From Q and one gets (see 

{a, b} = 4(v9"e • e + V""^^ ■ = 4(ci(x)a + {h{x) ■ Ob) (5) 

for some ci and li. Setting C = and ^ = Vip gives ci = — ^"^^^ and 
li ■ Vip = — ci. By taking the terms in which are of second order in ^ it 
follows that — cil — 2/1 (V(/?)*)^ • ^ = 0, which implies 



^" - cii - h{v^Y - v^ihY = 0. 

WW 



Applying this matrix to V(/? gives li = 



Now Vm = — 2m/i, and 

{ma, m6} = Vg(ma) • Vximb) — Vx{ma) ■ V^{mb) 
= m^{a, 6} + mb\7^a ■ \7x1ti — maSJxfn ■ 
= m^{{a, b} - Abh ■ ^ + Aah ■ Vip) = 0. 



□ 



We now specialize to the weight ^{x) = log and we will compute an 
explicit change of coordinates near T, which makes H^p into a d operator 
near S. In fact, the change of coordinates will consist of finding a codimen- 
sion 2 manifold in a neighborhood of S which is transversal to the flows of 
Hjna and H^f,- We then need to check that the flows originating from this 
manifold cover a full neighborhood of S in the cotangent space. 

Consider a truncated cone Vq = {x G R" : x„ > c|x|,c < |x| < c~^} for 
small c > 0. In suitable coordinates, Vq will contain Q. We will work in a 
fixed neighborhood 

V = V{6) = {xGR": x„> {c-d)\x\,c-6 < \x\ < c'^ + 6} 

with 5 > small. Now S is given by G T*V: x ■ ^ = 0, |x||^| = 1}, 

and we define a neighborhood 

U = U{6) ={{x,0 €T*V: \x-^\<5,l-6 < \x\\^\ < 1 + 5}. 

The first two new coordinates will be 

yi = x-^, y2 = \x\\^\, 

so S is given by {yi = 0, 2/2 = !}• 
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Wc have m = |xp and ma = \x\'^\(,\'^ — 1, mb = 2x ■ and the Hamilton 
vector fields are 

Hma = 2(|x|2e • V, - lepX • V^), Hrnb = 2(x • V, - ^ ' V^). 

These are smooth in T*R", and it is possible to compute the flows explicitly. 
If X and ^ are linearly independent, then the flows starting from (x,^) are 
given by 

Omais,ix,0) = (e2^^1x|((cos 2^y^ - yls)x + {sin 2^ yl - yls)Jx), 



e-2?^iiei((sin 2^yl - yfsp^ + (cos 2^y2 _ ^2,)^)^^ 

embit, {x,0) = {e^'x,e-^'0, 
where z = |||, and for z in the oriented 2-plane T = span(a;, ^), Jz is defined 

by 

Jz = —{z • ^)x + (z • x)^ 

where ^ = ||_||^^|^| ■ Thus, Jz is the unique vector in T for which \Jz\ = \z\ 
and (z, Jz) is a positive orthogonal basis of T. 

The leaf T^^^ through (x,^) G T*V is the set of all points in T*V which 
can be reached from (x,^) by the flows of Hj^a and Hj^b- 

Lemma 3.2. If (x,^) G U then F^^^ C U and 

^x,^ = {{z,v) ^T*V: z-i] = yi, \z\\ri\ = y2, span(2;, ??) = span(x, 0}- 

Proof. Denote the set on the right by L^^^. It is easy to see that T^^^ is 
contained in L^^^ by checking that z ■ r/, |2:||77|, and span(z,77) are constant 
along the flows. This also implies F^. ,c C U. 

For the converse, note that since x &V the plane T = span(x, ^) does not 
lie in {x„ = 0}, and we may deflne 

/ ProjT en 

w = w[x,^) - 



Iproj^^Cnl 

C = C{x, = yi^ + \Jyl-ylJ'u^- 

It follows that {w,C,) G Lx,e,- We claim that (x,^) = Omb{t,dmais,{w,Q)), 
provided that 

2vy2 - vi 

i = - log |x| - yis, 

where we deflne Z(e*") = a for — tt < a < tt. In fact we have 

(cos 2\^ y\ — y\s)w + (sin 2\^y\ — y\s)Jw = (x • ■u;)i/; + (x • Jw)Jw = x 

and so Oma{s, {w,C)) = (e^^^^x, •) and 9mbit,0ma{s,iw,C))) = {x, ■). The 
claim follows since any two points in L^^^ whose first components are the 
same must be identical. The construction also guarantees that the flow from 
{w,C) to (x,0 stays in T*V. 
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Now, if {z,r]) E L^^^, then w{z,r]) = w{x,£,) and Ci^TV) = Ci^^O- This 
shows that one may reach both (z, rf) and (x, ^) by flows starting from the 
same point {w^Q), which impUes [z,rj) S r^.^^. □ 

Since the leaves are given by the points where yi, y2-, and the plane 
span(x,^) are constant, we take the next new variable to be this oriented 
plane. More precisely, we take (ys, . . . ,y2n-2) to be local coordinates cor- 
responding to span(x,^) on the Grassmannian G{2,n), which consists of 
oriented two-planes in R"^. There is a single chart which achieves this, due 
to the fact that span(a;,^) does not lie in {xn = 0}: one may apply an ori- 
ented version of Pliicker coordinates, or if n = 3 it is sufficient to use the 
identification G{2, 3) = and stereographic projection. 

The final coordinates will be the flow variables y2n-i = s and y2n = t, 
where s and t were given in the proof of Lemma 13.21 The codimension 2 
manifold transversal to the flows will then be {?/2n-i = y2n = 0}. We have 
arrived at the desired change of coordinates. 

Lemma 3.3. The map ^ : {x,$^) y deflned in the discussion above is 
smooth and injective on U{5) for 5 small, and it is a diffeomorphism onto 
its image in R^". In the new coordinates Hmp becomes ^ + 

Proof. If <I>(x,^) = where (x, ^), (x', G U, then Lemma [3.21 im- 

plies that (x',^') is on the leaf through (x,^). Also, the vectors w,C and 
the coordinates s,t in Lemma 13.21 are the same whether they are computed 
from (x,^) or (x',^')- I* follows that (x,^ = 9mb{t,Gma{s, (w,C))) = {x',^,'), 
which shows that $ is injective. Since $ is smooth, it is a diffeomorphism 
onto its image. Further, Hmp becomes ^ -|- because the flows of Hma 
and Hmb commute. □ 



Note that Hmp = mHp on E, so the following result states in particular 
that one may solve HpU = / on S. 

Proposition 3.4. If /(x,^) G C°°{U{26)) for 6 small, then the equation 

HmpU = mf in U{5) 
has a solution u G C°°{U{6)) satisfying for all /c > 

|5°afu(x,0l < Ck\\f\\w^,oo^ui2S)), G U{5) 

whenever |a| -|- \P\ < k. 

Proof. Writing u = u o $, it is enough to solve 

+ idy,Ju = ixmf) o <^>-^ on 1>(t/(25)) 

where x £ C^{U{25)) and x = 1 ii^ar U{5). This may be solved using 
the Cauchy transform {dy^^_-^ + ^t^s/2n)~^' and the norm estimates are an 
immediate consequence. □ 
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4. CARLEMAN ESTIMATES 

In this section we will first recall the Carleman estimates for Hwy in 
and we note that the shifted estimate is valid when W e W^'^ D L'^{^}; C"'). 
We then use the Carleman estimate to solve an equation involving a conju- 
gated version of Hwy- 

The Carleman estimate in ^ is an estimate for the magnetic Schrodinger 
equation proved for limiting Carleman weights. Denote by (7 C R" an open 
set such that C ri, and recall that is a limiting Carleman weight in 
if (jnj and Q are satisfied in T*Q. Later we will restrict ourselves to the 
particular logarithmic weight. Below we write dil.± = {x G d^l: ^di,ip{x) > 
0}, and A<BiiA<CB where C > is a constant independent of h and 
£. Recall that || • || = || • ||l2(q). 

Proposition 4.1. Let he a limiting Carleman weight on Q. Suppose 
W G W^i'" n L'^{Q; and V G Then for ue H^n H^{Q) and h 

small, we have the Carleman estimate 

— h{d,^ipeh d,^u\eh d^u)QQ_ + ||eftti|| + ||eft/iVti|| 

<h\ehHy/yv}\ ^h{dyLpe~hdyu\e~hdyU)Q^^. (6) 

p __ 

Proof. By introducing v = e^u it follows that (joj) is equivalent to the a 
priori estimate 

- h{du(pduv\duv)gn_ + + ||/iVt;|p 

< h^We'hHwye'hvf + h{du^dyv\dyv)aQ,^ (7) 

for the conjugated operator ehHy^ye h . To prove ^ the idea is to work 
with the convexified weight 

if{x) = ip{x) + 

ip_ !£. . 

and first obtain the estimate for Hw,ve~ ^ ■ This operator can be split 
into a self-adjoint term P, a skew-adjoint term iQ, and a remainder term R, 

ehHwye'h = h~'^{P + iQ + R), 

where 

P = -h'^A- {V(pf, 

Q = hD oV(p + V(po hD, 

R = 2hW ■ {hD + iVip) + h^{W^ + D -W + V). 

We then have from ,5,, Eq. (2.12)] the estimate 

^{\\v\\^ + \\hVvf) <\\{P + iQ)v\\'' + h\d,^d,v\duv)en (8) 
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for some sufficiently small e > 0. By the assumptions on W and V, it follows 
from that 

\\Rvf < h\\\W\\l^^^^ + \\W^ + D.W + V\\l^^^^){\\vf + WhVvf) 

<h\\\vf + \\hVvf). 

Hence by choosing e > sufficiently small, the term can be absorbed 

in the left hand side of (jSJ and consequently 

1,2 

— {\\vf + \\hVvf) < h%^Hwye'^vf + h^{d^ipd^v\d^v)9n, 
from which ((Jj) can be derived as in [S]. □ 

Next we show that the Carleman estimate © can be shifted to a lower 
Sobolev index. For this we will use the semiclassical Sobolev spaces H^^^^RT') 
with norm ||/|U.^,(Rn) = UhDy fh^n-), where (0 = (1 + lei')'/'- 

Proposition 4.2. Let 99 be a limiting Carleman weight, and suppose that 
W E T^i'" n L°°(0; C") and V G Then for v G C^i^) we have the 

estimate 

II^IIl2(R") < Hw,ve-huf^^^,^^^y (9) 

Proof. Using the notation in Proposition 14.11 we have from (Sj Proposition 
2.4] the estimate 



/i2 



2 



for V G C^{Ct), where Cl is open and C f2 and C (7. Suppose we have 
extended W, V to R". The semiclassical counterparts of Q and ^ are 

II«^IIhi^j(R") ^ l|a||vKi'"nL°°(R")ll^ll_ffi^j(R")! 
and the corresponding dual estimates yield 

II^^IIh~1(R") ^ ^II^«IIl2(R")- 

scl ^ 

Hence for e > sufficiently small we have 
1,2 

''■ II ||2 ^ t4|| £ rr — ^/T 7-.\ Il2 

— Mh^jr'-)^'^ We'^Hwye h {hD}v\\ ^-.^1 ^^„y 

from which Q can be derived as in [S]. □ 

As a consequence of the estimate © and the Hahn-Banach theorem, we 
have the following solvability result. 

Proposition 4.3. Let (phea limiting Carleman weight, and let W G W^-'^'^H 
L°^(S7; C") and V G If h is sufficiently smah, then for any / G L'^{n) 

the equation 

eh Hw,ve~hu = f (10) 
has a solution u G /^^(n) with 

ll^i||L2(n) + II^-Du||l2(q) < /l||/||L2(n). (11) 
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5. PSEUDODIFFERENTIAL CONJUGATION 

In this section we will prove Proposition 14.31 in the case of uniformly 
continuous coefficients. This will follow by extending the pseudodifferen- 
tial conjugation technique introduced in ^H] (see also ^1]) to logarithmic 
Carleman weights 

ip{x) = log\x — Xo\, Xq ^ ^l. (12) 

The main result is the following. 

Proposition 5.1. Let (p{x) be as in dJ, and assume that W G (7(17; C") 
and V • VF, y G L"(17; C). If /i is sufficiently smah, then for any / G L'^{Q) 
the equation (fTUI) has a solution u G H^{Q) satisfying (fTTl) . 

To prove this, we start by choosing 17 to be an open set containing with 
xq ^ 17, and we take coordinates so that xq = and 17 C {x„ > c|x|,c < 
|x| < c-^} for c small. Note that Hwy = -A + 21^-L» + F with V G L"(17). 
We ignore the zero order term for the moment and consider the conjugated 
operator 

eh {{hDf + 2hW ■ hD)e~h = p + hQ 

where P = {hD + iV^pf and Q = 2W ■ {hD + iV(p). 

If h is sufficiently small. Proposition 14.31 (in the case W = V = 0) implies 
that we may solve Pu = / in 17 for any / G L^(17), and the solution operator 
P~^ : / I— > is a linear map which satisfies 

WP-'fWmn) + WhDP-'fW^,^^) < h-^f\\^,^^y (13) 

We would like to solve {P + hQ)u = / in the same way. If ||W^||j;^oo(q) is 
small then P + hQ is a small perturbation of P, and a solution is obtained 
from the Neumann series. If W is large a different method is needed. In 
this case we will find order pseudodifferential operators which conjugate 
the operator into a small perturbation of P. We use standard classes of 
semiclassical symbols and operators. 

Definition. If < o" < 1/2 and m G R, let S*™ be the space of all functions 
a(x,^; h) where x,^ G R" and h G (0, /iq], /iq < Ij such that a is smooth in 
X and ^ and 

\d^d^^a{x,^;h)\ < C7„^/i-'^l-+/5|(l + |Cp)™/2 
for all Q, /?. If a G we define an operator A = Op/^(a) by 

Af{x) = (2^)-" 11^^^^ e'^^-yy^a 

Note that we define the operators using semiclassical Weyl quantization. 
Operators in are bounded on with norm uniformly bounded in h, and 
the composition of two operators is again an operator in the same class. 
We assume familiarity with semiclassical calculus in what follows, for more 
details see jl] and also ^3 where the result of this section was proved in the 
semiclassical setup for linear Carleman weights. 



)f{y)d^dy. 
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To deal with the nonsmooth coefficients we extend W to a vector field in 
Cc{^', C^) and V by zero, and we consider the mollifications 

W^ = W*X5, V^ = V*xs (14) 

where 5 = h'^ and < o" < 1/2. Here x&{^) = ^~^''x{x/^) is the usual 
mollifier with x G C^^iW), < x < 1, and j xdx = I. We write W'' = 
W — and V'' = V — and note the estimates 

\\d"W^\\Loo, ||5"y»||Ln = ©(/i-'^l"'), (15) 
\\W^\\l-°,\\V^\\l" =o{l) (16) 

as /i — > 0. 

Write = 2W^ ■ (^ + iVtp) - hD ■ W^, so that Q« = Op^iq^) = 2W^ ■ 
(hD+i'Vif). The following lemma shows that one may use pseudodifferential 
operators to conjugate away the main part of the first order term Q. 

Lemma 5.2. There exist c, c, r G S'^ so that 

(P + hQ^)C = CP + h^'^'^R in n. 

Both C and C are elliptic, in the sense that c and c are nonvanishing for 
small h. 

Proof. Suppose c is any symbol in 5^ which is equal to 1 outside a large ball 
in (x,^). Then for some tq G 5^, 

(P + hQ^)C = CP + hOpJ-HpC + q^c) + h^-^''Ro in (17) 

We will choose c so that jHpC + q^c is of lower order. Trying c = e**^ and 
using the notation in Section |S1 and Proposition \\\A\ we choose to be a 
solution of 

Hmp(l) = -mi in [/ = \J{6). (18) 
Since x and ^ are bounded in f/, (fT5|) gives 

Let be a smooth cutoff supported in C/((5) which is equal to 1 on 

U{5/2). Then c is chosen to be 

c = e'^^. 

It follows that c G 5^ and c is nonvanishing. 
Since Hmp = mHp + pHm, we get by (|18() 

m{^HpC + gt'c) = {{HmpX)(t> + xHmp(t> - pHm{x4>) + mq^)c 
= {{HmpX)4> + (1 - X)mq^ - pHm{x4>))c. 
Let ip G C^(il) and ip = 1 near Then ^{\HpC + qh) = Ip where 

I = + J-(H^^x)<t>--H„,{X(t>))c 

\ p mp m J 

and / G h^^S^. Then (dZJ gives 

{P + hQ^)C = CP + h^-^''R inn 
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where c = c + hi is of order and elhptic, and r G S^. □ 

We may now prove the main resuh. The proof involves a number of 
cutoffs to restrict the functions to the set O where if is smooth. 

Proof, (of Proposition 15.11) Using the notation in this section, we want to 
solve 

{P + hQ^ + hQ^ + h'^vH h^V^)u = f inn (19) 

where = 2W^ ■ {hD + iVip). Let V'i {j > 1) be C^{n) functions with 
V'l = 1 near and Tpj+i = 1 near supp{tpj). We try a solution 

u = Cip2P'^^3C-^Tp^v (20) 

for some v £ L^(R"). Here is the inverse of C on L^(R"), which exists 
for small h. Inserting (|20j) in (|19|) and using Lemma 15.21 gives 

{M + Ri)v = h^f inn (21) 

where M = C'PV2^'" Vsf?" V4 and HV'i-RiIIlz^l^ = o(l) as /i ^ 0. Here 
we have used (UJ), (UHI), (CSl), and the identity djC = Cdj + Op;,(9a;.c). 
We may write 

M = tlj2 + [C, V'alC'- V4 + C[P, ^^2]^- VsC" V4. 

Extending / by zero and modifying (|2H) slightly, we obtain an equation in 
R" of the form 

{I + Ml + ^iRi)v = f inR" (22) 

where Mi = jpi[C,jp2]C-^iJ4 + ipiC[P,ip2]P~^^3C~^iJ4- By the pseudolocal 
property ||Mi||£2^j^2 = o(l) as h ^ 0, so we obtain a solution f G L^(R"') 
to (|22|) by Neumann series. Since ipj = 1 on Q the function v will satisfy 
H21|) near $7. It follows that ((20} solves H19() . and one obtains the desired 
norm estimates from H13() . □ 



6. Construction of solutions 

We now give the construction of special solutions to Hwyu = in 17. 
This proceeds as in 0, except that an additional smoothing argument is 
required. The solutions will have the form u = e~hv where f is a WKB 

ip ip 

solution for the conjugated operator eh Hw,ve~''^ ■ Thus, u is of the form 

u = e~h{a + r) (23) 

where p = ip + and tp is a real valued phase function, a is an amplitude, 
and r is a correction term. By Lemma 12.21 we may assume that V • = 0, 
and then 

h^Hwy = {hDf + 2hW ■ hD + h^V 
with V € -L"", and inserting (|23j) into h'^Hwyu = gives the equation 

{{HD + iVpf + 2hW • {hD + iVp) + h'^V){a + r) = 0. 
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Collecting like powers of h, u will be a solution provided that one has in 17 
{Vpf = 0, (24) 

{Vp-D + Vp-W + ^/^p)a = Q, (25) 

ehh?Hwye~hr = -h^Hwya. (26) 

We fix </?(x) defined by (|12j) . Now (|24|) is an eikonal equation for ^, which 
reads 

{V^l:f = {V^f, V(/9-VV' = 0. 
As shown in ^9-^, a solution is given by 

^{x) = dist5„-i ( ,^~^V u;) (27) 

\\X — Xq\ J 

where uj G is chosen so that iIj is smooth in (l. 

We choose coordinates so that xq = and $7 C {x„ > 0}, and we take 
u = ei. We also write x = {xi,r9) with r > and G 5""^^, and consider 
the change of variables ^' : x i-^ (z, ^) where z = xi + ir is a complex 
variable. Writing / = / o ^f"-*^, we get p = log z, (Vp) = ^(ei + icr) where 
er = (0,9), (ApY = and Vp • V becomes fO,-. 

We see that (|25|) is a 9-equation in the new coordinates, and the solution a 
will have in general the same regularity as W. On the other hand, in (|26j) one 
needs two derivatives of a on the right. To deal with this for W £ C^{^1; C") 
we use the same approximation procedure as in (|14|) . now with a > small 
and W extended as a vector field in C^{il; C"). The solution will be taken 
of the form u = e~h(a + r), where p is as above and one has in 

{Vp- D + Vp-W^ + —Ap)a = 0, (28) 
2i 

ehh'^Hwye-Tir = -2ih{Vp ■ W^)a - h'^Hwya. (29) 
To solve (|28|) we take a = a o so a must satisfy 

2 2[z — z)J 

Here = W"^ o We try d = {z — z)^~e**, and we get the equation 

= -^(ei + ier) ■ in ^(O). (30) 



It is now easy to solve (j2^ 
Lemma 6.1. The equation (|28)) has a solution a E C°°(r2) which satisfies 

Proof. We first take x ^ with x = 1 near Q, and consider 

5^$ = -^X(ei + ier) • W^" for z G C 

with X = X ° This has the explicit solution 

$(z, 0) = — 3- [ -x{z - w, 9){ei + ier) ■ W^{z - w, 9) dw A dw. (31) 
Ani Jq w 
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Then a = {{z — z)~2~e**) o will be a solution of (|28]) in 17 with the given 
norm bounds. □ 

It remains to solve (|29|) . The norm of the right hand side of (|29j) is 
0{h^~^'^^) when a is small enough, using the Holder continuity of W. It 
follows from Pr op osition 15 . II that when h is small there is an H^{n) solution 
r with ||r||^2(f,) = 0(/i-^), ||Vr||^2(f,) = 0(/i-i+-^). 

We collect the results obtained in this argument. 

Proposition 6.2. Let W £ C^(J7;C"), V • 1^ = 0, and F G L"(fi). Let 
ip, ip be defined by (|12|) and (|27|) respectively. Then for h small there is an 
H^(Q) solution u = e~h^'^~^'^'^\a + r) of the equation H]yyu = in 0, where 
a is given in Lemma l6.ll Also, one has the norm estimates 

||a°a||i^(f,) = 0(/i-H), 

Mmn) + \\hVr\\L2^n) = Oih^') 

where a > is small. 

7. Determining a magnetic field 

In this section we prove Theorem 11.11 following the arguments given in [2]. 
Assume that Wj G C%n; C"), V • Wj = 0, Vj G and that 

Nw^,vJ\f = Nw,,vJ\p for all f G H^/\dn). (32) 

Also assume that Wj are extended as vector fields in C^{Q; C"). 

We start by choosing ui and U2 to be solutions of H]y-^y^ui = and 
Hyy^ y^U2 = 0, of the form 

where ^{x) and tpi^x) are defined by (|T^ and respectively, Uj satisfies 
(V((/p + #) • ^ + V((^ + iV) • Wl + ^A((^ + = 0, 

(V(-(^ + + V(-(^ + • + ^A(-(^ + #))a2 = 0, 

and ||9°aj||ioo(f^) = 0(/i"'"l"l) for cj > small, and finally ||rj|| = 0(/i'^^), 
llVrjII = 0(/i-i+'^^) as /i ^ 0. 

Take tt2 to be the solution to ffvi^2,V2^2 = with U2 = ui on and let 
u = Ul — U2- We use Lemma 12.31 for u and ii2 to get 

{Hw2y2u\u2) = -{duu\u2)dn- (33) 
We examine the right hand side of (|33|) . 
Lemma 7.1. One has /i(9jyii|n2)ao — > as /i — > 0. 
Proof. Because of (|32]) we have d^u = i{{W2 — Wi) • i')ui on F, and so 

h{duu\u2)dn = ih{{{W2 - Wi) ■ i^)ui\u2)p + Kduu\u2)g^^p (34) 
The first term satisfies 

h\(,{{W2 - Wi) ■ i^)'Ui|'U2)^| < h\\ai + ri||i2(gQ)||a2 + r2||i2(gQ). 
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We have ||aj||i2(gn) < llL'-(an) = 0(1) and \\rj\\L2(^g^) < \\rj\\Hi/2+s(fi) = 
0{h~^/'^~^"'^~^) for any 6 > 0, which gives 

\\aj+r,h'^dn)=Oih-'/^+') (35) 

for some new 6 > 0. This shows that the first term in (|34|) vanishes as 

For the other term we compute 

To estimate the normal derivative we use the Carleman estimate of Propo- 
sition for the weight —ip, in the form 

Vh\\y^d,yipe~h d^u\\i2(^g^^^i + ||e~ftii|| + ||e~ft/iVu|| 

< h\\e~hHw2,V2^\\ + ^\W-9uVe~hduu\\L2(^Qf^_y 
The estimate apphes for this u since u £ H/\{^}) and u\gQ = 0. Since 

e~hdvU = ie^^{(W2 — Wi) ■ v){ai + ri) on F, and since dO.- C F, it follows 
from that 

<Vh\\e-^Hw2y2u\\+0{h-^'^+'). 
Using that ui and U2 are solutions and V • Wj = 0, we obtain 

Hw2y2U = Hw2,V2Ui = 2{W2 - Wi) ■ Dui + [W^ - Wf + V2- Vi)ui. 
The explicit form for ui implies 

-^1^2,^2^^ = e^^'P^'^\2ih-\Wi - W2) ■ V(v9 + ii))]{ai + n) 

+ 2{W2 - Wi) ■ Dim + ri) + (W| - Wf + V2- Vi)(ai + n)) (36) 

which shows that ||e~ -ffyi/jya^ll — 0{h~ ), since the L"-2 norm of ai + ri 
can be estimated by the norm. Collecting these estimates gives 

h{d^u\u2)9^^p = 

as h ^ 0, which concludes the proof. □ 
From (|36j) we obtain 

h{Hw2,V2'^\u2) = / 2i{Wi - W2) -Viip + iij)aia2dx + o{l). 
Jn 

Then (|33|) and Lemma l7. 1 1 implv 

/ V(v9 + ii/j) ■ (Wi - W2)a dx = (37) 
Jn 

where a = lim/i^o '2ifi2- We will now switch to the complex notation as in 
Section El 

Choose coordinates so that xq = and uj = ei, and write x = {xi,x') 
where x' = r9, r > 0, and 6 G 5"~^. Let Pg be the two-dimensional plane 
consisting of points {xi,r9) for 6 fixed, and write = 17 n Pg. We also use 
the complex variable z = xi + ir, which identifies Pq with C. 
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In the coordinates {z,0) = ^(x), Lemma l6 . 1 1 shows that a in (|37() is given 
by a o = (2 — ^)^~"'e**, where $ satisfies 

d,^ = -^{ei + ier)-{Wi-W2) m^{Q). (38) 

We have written Wj = Wjo^i/^'^ . Note that since ^ is of the form (|ST|) . Wj G 
imphes that ^> G and further ^>( • , 61) E (7^+^(0). The latter follows 
since 9^$ G and G d^d^^C^, and the singular integral operator 
dzd^^ is bounded on Holder spaces. 
Now ((H7j) becomes 

/ (/ -{ei + ier)-iWi-W2)e''^dzAdz)de = 0. (39) 

We need a slightly more general result. From the transport equation for 
oi we see that in the definition of u\ we may replace a\ by a\g\^ where 
V((/9 + i-^) ■ Vgi = and gi € W'^'°^{Q) is independent of h. It follows that 
(|39|) holds with e** replaced by e**gi. Choosing = zg{z)g{6) with Sfi^ = 
and g smooth, and by varying g, we see that for almost every 9 G S""^ one 
has 

/ (ei + ier) ■ (Wi - W2)e''^giz) dzAdz = 0. (40) 

The argument that ()40() implies dWi = (iW2 now proceeds as in However, 
this was written in 5 quite briefly, especially in the case where has 
nontrivial topology. Therefore, we will give a rather detailed argument for 
determining the coefficients. 

It will be enough to restrict to 3-dimensional subspaces L of R" which 
consist of the points {xi,r6) where 6 G S""^ varies on a fixed two-plane. 
The fact that (HOI) holds in each L wiU imply dWi = dW2 in 0, see 0. Thus 
we may assume n = 3 in what follows, and 6 ^ may be identified with 
an angle, using Ze*° = a when < a < tt. 

We begin by looking at what kind of topology VIq can have when 9 varies. 
For this we use the map : dVL — > R, x 1-^ '^W\' "^^^ main point is that 
possibly after a small change of coordinates, is a Morse function. 

Lemma 7.2. For almost every choice {u!,xo) G TS"^, is a Morse function 
and there are finitely many critical values ^0 < ^1 < • • • < ^iv of 0. Also, 
e{dQ) = [60, On], and if ^ G {Oj,ej+i) then 

(a) do, and Pg intersect transversally, or equivalently for any x G dQCiPg 
the vector I'ix) is not orthogonal to Pg, 

(b) Qg is a bounded open set with smooth boundary in Pg, and dilg = 

dnnPg, 

(c) there is a diffeomorphism T : {9j,6j^i) x dOg Q~^{{9j,9j^i)), 
and if Tt{x) = T(t,x) then Tt is a diffeomorphism of di^g and d^lt 
with Tg = IddQg . 

Proof. We see that x G 90 is a degenerate critical point of if and only if 
N{x) ■ Lu = N(x) ■ {x — xq) = and K(x) = 0, where is the Gauss map 
and K is the Gaussian curvature of dQ. Thus, is Morse provided that 
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(w, xq) is not a critical value of the map F : ^7(517) x R ^ TS'^ 
F{{x, Lo), t) = {uj, {x ■ N{x))N{x) + t{uj X N{x))), 



where Q{dQ) C T{d^}) is the unit sphere bundle. Sard's lemma shows that 
is Morse for almost every (a;, xq). Then Q has finitely many critical values, 
and if 6 is not a critical value then dil and Pq intersect transversally. 

The other condition in (a) is just another way of stating the transversality, 
and this latter condition implies that OQq = di^DPQ. Then, OQq is a smooth 
1-manifold and 0,q has smooth boundary in Pq. Part (c) follows from the 
fact that there is a gradient-like vector field for Q 11 , and we obtain the 
diffeomorphism T by considering flows of this vector field which originate 
from O^Iq. □ 

Note that d^le can have many components, but if 9 stays away from the 
critical values then the number of components stays fixed and the compo- 
nents vary smoothly with 6. Unless stated otherwise, we will assume that 
the coordinates are chosen as in Lemma 17.21 and 9 is not a critical value. 

Lemma 7.3. dWi = dW2 in 

Proof. Our starting point is ()4U|1. Because of the factor e** this may be 
considered as a nonlinear Radon transform, evaluated at the plane Pg. Most 
of the work below is to show that ()4U|) remains true with e** replaced by 1, 
which corresponds to the usual Radon transform. 

The proof is in several steps. For the complex analysis terminology see 



Step 1. Using the equation (|H5|) and integrating by parts, we obtain from 
H4U|1 the orthogonality condition 



for any holomorphic g G W'^'°°{Qe). Since Pg \ has finitely many compo- 
nents, an approximation argument (using a version of Mergelyan's theorem) 
implies that (|41|) holds in fact for any holomorphic g E C{VtQ). 

Step 2. The condition H41|) is equivalent with saying that e^^lgj^^ is the 
boundary value of a holomorphic function F G C{Vl0). Indeed, write / = e**, 
and let F be the Cauchy integral of f\dng- From (|^T|) we see that F = 
outside Og, and the Plemelj-Sokhotski formula then implies that F G C{Q,e) 
with Fla^e = /• 

Step 3. We would like to show that F is nonvanishing and has a holomorphic 
logarithm in 0,g. We first claim that if 7 is a closed curve in dQg, then 



To show this we write fs = e*'^* for < s < 1. We obtain that / o 7 is 
homotopic to /o o 7 = {1} in C \ {0}, and the claim follows. 

Step 4- We can use H42() and the argument principle to conclude that F is 
nonvanishing, also when fig is not simply connected. Let 7 be an oriented 




(41) 




(42) 
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parametrization of OQq as a sum of simple closed curves. It follows that 
Ind-y(Q;) = Jg^^ dz = for all a ^ ftg. Also, F 7^ on Sf^g, so there are 
only finitely many zeros in ^Ig. One may now use the residue theorem and 
argue in the usual way that ^ fp^^ ^ dz = is the number of zeros of F in 

Step 5. Next we will show that F has a holomorphic logarithm in Qg. Again, 
this would be immediate in a simply connected domain. In the general case, 
F has a holomorphic logarithm provided that 

^dz = (43) 

for any closed curve 7 in ^Ig. To show ()43() for given 7, we take Cj to be 
the finitely many components of P51 \ Jig, and for each j we let 7^ be a cycle 
corresponding to the oriented boundary of Cj. Then 7^ is contained in d^g, 
the index of jj on Cj is 1, and the index of 7^ on any other component 
is 0. Thus, using the 7^, we can construct a cycle 7 contained in dQg so 
that the index of 7 is equal to the index of 7 at each point outside of ^Ig. It 
follows that 

F' f F' [I 

dz = — dz = - dz = 



using (|42|) . This shows (|43)) . and we obtain a holomorphic logarithm G G 
C{Qg) with F = by fixing a point zq in each component of Qg, and by 
taking in this component 

Giz)= r^dw + co (44) 

where the integral is over any curve connecting zq to z, and e'^ = F{zo). 

Step 6. Since e*^"** = 1 on dQ.g, we get Gla^g = i<I> + w where v is constant 
on each component of dVtg. In fact, v is equal to 27rim for some m G Z on 
each of these components. This shows that 



Gdz= i^dz. (45) 

Step 7. We now return to (|1T|) and take 5 = ^ G C{i}g). Then (|45|) gives 

^dz = 0. 



'ang 

Integrating by parts and using yields 



/ 



(ei + ier) ■ {Wi - W2) dzAdz = 0. (46) 



Step 8. We need to show that is valid with er replaced by — e^. If Wj 
were real valued this would follow just by taking complex conjugates. If the 
Wj are complex, we can go back to the beginning of this section and repeat 
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the construction of solutions, with ip replaced by —ip. In this way, instead 
of (|57j) we arrive at 

/ V{ip + iip) ■ {Wi - W2)adx = 
where a o '^^^ = — z)^~"e**, and <I> satisfies 

d,^ = \{ei+ ier) ■ {^i - 1^2) in ^{^). 

Going through steps 1 to 8 above, we obtain with Wj replaced by Wj, 
and taking conjugates gives 

/ (ei - ier) ■ {Wi - W2) dzAdz = 0. (47) 

Step 9. If ^ G Pe then ^ is a linear combination of ei and e^, and we get 
from (011) and (jUj) that 

[ ^-{Wi- W2) dzAdz = 0. 

JUg 

Returning to the x coordinates, this gives 

[ ^-{Wi- W2) dS = (48) 

J(xo+P)nn 

for all two-planes P passing through ei, and all £ P. 

Step 10. The left hand side of (|48|) is related to the Radon transform of an 
expression involving curl (Vl^i — W2). Varying xq in a small neighborhood of 
and varying a; in a conic neighborhood of ei, one obtains curl {Wi — W2) = 
from the arguments in Lemma 5.2]. □ 

After showing that dWi = dW2, the final step is to show that Vi = V2. 
This also follows along the lines of [Sj, but we need to give more details to 
account for the nontrivial topology of the 

Lemma 7.4. Vi = V2 in Q. 

Proof. The proof is again in several steps. 

Step 1. Since curl(iyi — W2) = and is simply connected, we have 
W2-Wi = Vp where p G (7^+^(11). 

Step 2. The assumption that N-\j^^y^f\p = N-\ji/^y^f\p for all / G H^/'^{d^), 
together with boundary determination results for the magnetic Schrodinger 
equation P, implies that (M^i)tan = (W^2)tan on F. Note that the results in 
^ remain valid for V G L"(r2). It follows that the tangential derivatives of 
p vanish on F, which shows that p is constant on each component of F. 

Step 3. We use again the complex notation, with xq = and uj = ei. 
Inserting W2 — Wi = Vp in (|l-i8|) . we see that ^ = p will be a solution 
of ()38() ii p = p o Repeating the arguments in the proof of Lemma 

17.31 shows that e*^|gQg is the boundary value of a nonvanishing holomorphic 
function, and taking a holomorphic logarithm gives ip = Gg + vg on dO,g, 
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where Gq G C{^q) is holomorphic and vg is equal to some Inim, m S Z, on 
each component of dflg. 

Step 4- We claim that for each component Cq of Qq, there is a point yo £ dCg 
with uq £ F. In fact, we may take j/o to be a point which minimizes \y\ among 
y £ Cg. Then yo £ dCg and (1 — t)yo ^ Cg for t > 0. It is enough to show 
that 2/0 ■ t^iuo) < 0, since then yo S F. But if one had yo ■ I'iyo) > 0, then 
the fact that for any c > there is a truncated cone 

{y G R": (y - yo) • z^(yo) < -c|y - yo|, |y - yo| <S}nng 

contained in Cg, would imply that (1 — t)yo G Cg for some t > 0. This is a 
contradiction. 

5iep 5. From Steps 3 and 4 we see that Cg £ C{Qg) is a holomorphic 
function in Qg, and for any component Cg of i}g there is an open set of dCg 
in which Cg is constant. This implies that Cg is constant on each component 
of r^e, hence also on each component of di^g. 

Step 6. We also need that Cg\Qng can be made to vary continuously with 
9, as long as 6 stays away from the critical values of Q. To do this we use 
Lemma 17.21 (c) which implies that d^g varies continuously with 6, and so 
does F in (|44|) . Using the regularity properties of <I> = p, and choosing the 
point zq in (|44|) to lie in dVtg with continuous dependence on 9, it can be 
checked that Cg\Q^g varies continuously with 6. 

Step 7. We now restrict to the set Q~^{{9j, 9j^i)), and claim that p is locally 
constant in this set (which implies that p is constant on the components of 
this set, since p is continuous). To prove this, we first recall that ip = Cg+vg 
on dQ.g, where p and Cg vary continuously with 9. Then also vg varies 
continuously with 9, and since -^vg is integer valued we see that vg is 
locally constant. 

If X £ Q^^{{9j,9j+i)) and ©(x) = 9, then x lies in some component Cg of 
By Step 4 there is some point y in Ce n F, and Step 2 shows that p is 
constant in B{y, r) n dQ for some r > 0. The same then applies to Cg near 
y. Since the components of dVlg vary continuously, and since Cg is constant 
on each such component, we see that Cg is constant near x when 9 varies. 
This shows the claim. 

Step 8. To show that p is locally constant on dVl, it remains to check that this 
is true near each x £ dVtr\Pg where is a critical value. If ^{x) is orthogonal 
to Pg, then X £ F and this follows by Step 2. If i'{x) is not orthogonal to 
Pg, then for some r > and j, both the sets B{x,r) n Q~^{{9j,9j+i)) and 
B{x,r) n Q^^{{9j^i,9j)) are nonempty and connected. Since p is constant 
on both these sets and continuous, p must be constant also near x. 

Step 9. The preceding step implies that p is constant on the components 
of do.. This shows that (Wi)tan = (W^2)tan on do.. Recalling that 30 is 
connected, we may substract a constant from p to obtain p\dQ, = 0. The 
assumption, gauge invariance, and Step 1 then imply 

^Wiy^F = ^W2,V2\f = ^W2-vpy2\F = ^Wxy2\F- 
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Consequently, we may assume Wi = W2 in the arguments in this section. 
Going through the proof of Lemma 17.11 and using the assumption Vj G 
L°°{fl), we may take hmits as /i — > in H33() to obtain 

/ {Vi-V2)adx = (49) 
Jn 

where a o vi/^i = (^z — z)^"". 

Step 10. We argue as in Lemma 17.31 and replace a in by where 
V((^ + iip) ■ Vgf = 0. Moving to the variables {z,9) and taking 5 to be a 
function of 9, and by varying g, we obtain 

[ {Vi-V2)i^'\z,e))dz Adz = (50) 

for almost every 6. Here we use the fact that — V2 is in L^(R"), so the 
restriction Vi — V2|pg is integrable on Pg for almost every 6. 

In the X coordinates, (jSU)) reads 

[ {Vi - V2) dS = 0. (51) 

J{xo+P)nn 

This is valid for almost every two-plane P passing through uj = ei. Varying 
xq and (jj slightly, the Radon transform arguments in give Vi = V2. □ 



8. Determining a convection term 

Finally, we prove Theorem ll.2l bv reducing the inverse problem for the con- 
vection equation to the corresponding problem for the magnetic Schrodinger 
equation. 

Proof, (of Theorem II. 2j) We observe that Hiw^^q^ = —A -|- 2Wj ■ V where 
Qj = Wj — V • Wj. Consequently, for / G H^/'^{d^) we have the relation 
NiWj,qj f = ^Wjf — {Wj ■ z/)/ between DN maps, and 

{Niwu,, - Niw,,,,)f\F = -iiWi - W2) ■ y)f\p 

using the assumption that Nw-^f = Ny^/^f on F. By this assumption 
also implies that Wi = W2 on F, and we obtain [Nmr^^q^ — A^iVK2,g2)/l p 

Now Theorem 11.11 shows that dWi = dW2 and qi = q2 va. Since 
is simply connected there is p € C^~^^{^) with W2 = Wi + Vp. Theorem 
11.11 also shows that (TVi)tan = (W^2)tan on dVL, so p is constant on dVL which 
was assumed to be connected. By substracting the constant we may assume 
that p\dQ, = 0. The condition qi = q2 then implies that 

-Ap + {2Wi +Vp)-Vp = Q in 

Since everything is real valued, the maximum principle shows that the only 
solution with zero boundary values is p = 0. It follows that Wi = W2. □ 
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